Dynamic mechanical response of polymer networks 

S. F. Edwards, H. Takano* and E. M. Terentjev 

Cavendish Laboratory, University of Cambridge 
Madingley Road, Cambridge CB3 OHE, U.K. 

February 1, 2008 



Abstract 

The dynamic-mechanical response of flexible polymer networks is studied 
in the framework of tube model, in the limit of small affine deformations, 
using the approach based on Rayleighian dissipation function. The dynamic 
complex modulus G*(u) is calculated from the analysis of a network strand 
relaxation to the new equilibrium conformation around the distorted prim- 
itive path. Chain equilibration is achieved via a sliding motion of polymer 
segments along the tube, eliminating the inhomogeneity of the polymer den- 
sity caused by the deformation. The characteristic relaxation time of this 
motion r e separates the low- frequency limit of the complex modulus from the 
high-frequency one, where the main role is played by chain entanglements, 
analogous to the rubber plateau in melts. The dependence of storage and 
loss moduli, G'(lo) and G"(to), on crosslink and entanglement densities gives 
an interpolation between polymer melts and crosslinked networks. We dis- 
cuss the experimental implications of the rather short relaxation time and 
the slow square-root variation of the moduli and the loss factor tan 5(lo) at 
higher frequencies. 

PACS numbers: 

83.80.Dr Elastomeric polymers, 

62.40.+i Anelasticity, internal friction, stress relaxation, 
83.50.Fc Linear viscoelasticity 



'Permanent address: Department of Physics, Faculty of Science and Technology, 
Keio University, Yokohama 223, Japan 



1 



1 Introduction 



Viscoelastic properties of dense polymer melts have been successfully ex- 
plained by the reptation model approach, where one examines the con- 
strained motion of a chain within an effective tube formed by entanglements 
with other polymers jl], [2| . In this model, the effect of entanglements is to 
restrict the lateral motion of a chain and is modeled, in the first approx- 
imation, by a dynamic constraint in the form of a tube surrounding the 
thermally fluctuating polymer. The dynamics of the polymer is then de- 
scribed by the reptation motion within such a tube, essentially a diffusion 
of segments along the random primitive path, with an effective number of 
chain segments constrained by the tube a new dynamic variable. Accord- 
ingly, the characteristic disentanglement time for a chain is that required 
for all its length to diffuse out of the initial tube. On the other hand, the 
classical theory of rubber elasticity, see for example Treloar's treatise ||] 
or 0, examines the response of a crosslinked polymer network to an affine 
deformation in the opposite limit of high dilution: Only the topologically 
quenched random crosslinks are affected by the strain, while the rest of the 
chain is assumed free to explore all its configurational space. Several pen- 
etrating attempts have been made over the years to improve this simple 
approach and address the issue of dense, entangled network ^. Although 
the disentanglement by reptation out of its tube is impossible in a perma- 
nently crosslinked network, the concept of a dynamic primitive path, around 
which the polymer segments are constrained to fluctuate, has been used to 
describe the static elastic properties of entangled networks. The sliplink 
model Q, where the entanglements are treated as discrete overdamped but 
nevertheless mobile constraints, is also a variant of the tube approach. 

The main success of reptation models in polymer melts has been the de- 
scription of dynamic properties and viscoelastic relaxation rates. Analogous 
approaches to crosslinked networks are not well developed at the moment 
and, in fact, even the basic dynamics of stress relaxation in rubbers is not 
described theoretically. Much progress has been achieved in studies of sec- 
ondary effects, such as the relaxation of free dangling ends attached to the 
network, leading to a slow power-law decay of stress [||, ||]. A significant 
contribution to understanding of the role of entanglements in rubber elas- 
ticity is made by computer simulation studies, e.g. [10]. However, on a basic 
level, one first needs to analytically describe the relaxation of an affinely de- 
formed ideal network. In this paper we investigate the dynamic mechanical 
response (the linear complex modulus) of rubbery networks in the dense 
limit of high entanglement. One should note that, unlike in a case of free 
chains, the entanglements of networks strands ought to be understood in a 
more general sense. In particular, the knots of increasing complexity, de- 
scribed by corresponding topological invariants, e.g. [11], are not the only 
way to restrain the lateral motion of network strands: Even a simple loop, 
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both ends of which are held by crosslinks, will contribute to the formation of 
an effective reptation tube. Therefore, the traditional notion of minimal arc 
length required for a polymer to be entangled, cf. J2|, should not be applied 
to networks - even relatively short strands could be highly constrained in 
a dense system. In such networks, one expects to find a storage modulus 
G' even at zero frequency of imposed deformation (the static modulus /xo, 
determined by the crosslink density in the classical rubber elasticity ||), 
while at a finite frequency u> the response should be increasingly controlled 
by chain entanglements in the similar way as the rubber plateau value of G' 
in a dense melt. 

Let us consider a polymer chain with segments of Kuhn step length b, of 
which both ends are crosslinked. The entanglement effect of other network 
strands is represented by a surrounding tube of diameter a. The central line 
of such a tube, regarded as a sequence of taut paths between entanglement 
points, is given a name of primitive path and may be considered as a random 
walk of step length a, with the end to end distance the same as that of the 
polymer chain. In the equilibrium state, the density of the polymer segments 
distributes uniformly along the primitive path. The linear viscoelasticity 
is then described by the response of chain segments to the instantaneous 
deformation. As a starting point, one may assume that the deformation 
of the primitive path is affine, in the sense that both the chain crosslinks 
and the entanglement points move in proportion to the macroscopic strain 
imposed on the body of rubber. Immediately after such a deformation, the 
density of polymer segments along the primitive path is inhomogeneous, 
because the deformed primitive path steps have different lengths depending 
on their directions and the chain cannot respond instantaneously to the 
change of its constraints. In time, the inhomogeneity of the polymer density 
relaxes through the sliding motion of the chain along the new primitive path. 

In a network, the reptation motion of the chain is completely prohib- 
ited by the crosslinking. Therefore, this sliding motion of chain segments is 
considered to be more important for the dynamic properties of crosslinked 
networks than for corresponding melts, where the main effect is the repta- 
tion diffusion out of the deformed tube. In this paper, we clarify how the 
slip motion of polymer segments along the primitive path affects the dy- 
namic modulus. The results are summarised in Fig. |l](b). We concentrate 
on relatively low frequencies, well below the glass region. One expects the 
linear static modulus hq = G'(u —* 0) to increase with the network crosslink 
density, but also with the density of entanglements. In the limit of higher 
frequency, r e uj 3> 1, with the characteristic time T e of diffusion along the 
primitive path, the storage modulus approaches the plateau value corre- 
sponding to the rubber plateau in a melt, where the role of entanglements 
is dominant. Our challenge is to develop a universal description that would 
also incorporate the limit of zero frequency, and the plateau modulus in 
a dense melt. 
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(a) (b) 



Figure 1: (a) The scheme of tube model and the relevant variables: the 
polymer conformation r(s,t) and the tube primitive path R(n,t), where 
the function n(s, t) specifies the position along the fixed-length path and s 
specifies the position along the chain. In a network, the chain is permanently 
crosslinked at both ends of the tube. (b) The "expected" behaviour of 
the dynamic storage modulus G'{uj): For an uncrosslinked melt the static 
modulus is zero, but reaches a rubber plateau at a finite frequency due 
to entanglements. On increasing concentration of permanent crosslinks the 
static modulus G'\ w= q = (jlq increases; the plateau modulus has contributions 
from both crosslinks and entanglements. At a much higher frequency all 
curves should converge to the glass plateau, which is not the subject of this 
paper. 



A classical theoretical description of polymer networks is based on the 
statistical- mechanical description of Gaussian chains with topologically quenched 
constraints of permanent crosslinks, in the simplest case, and entanglements 
on a more advanced level j|, [Tj. The computation of a Gibbs partition 
function Z = J is insufficient because of the quenched distribution of 
crosslinks and the Replica method has been developed as an extension of 
classical statistical mechanics, as one of the ways to evaluate average macro- 



scopic parameters, such as the free energy -F expt = (log Z) [[|, 12]. However, 
in spite of many successes of the Replica method, it remains unsuitable for 
finding the dynamic properties, such as the complex modulus G*(uS) - a tra- 
ditional parameter describing viscoelastic behaviour of polymers and rubbers 
[13] . One needs a dynamical, not static equilibrium solution. Here we utilise 
a method of solving the dynamical problem of rubber-elastic response by an 
extension of Boltzmann approach, based on Langevin and Fokker-Planck 
equations. One introduces the Rayleighian 1Z = L — F + ^AC to incorpo- 
rate the chain equilibrium properties via its Lagrangian L, the kinetic effects 
via the dissipative function F (not to be mixed with the free energy!) and 
all required constraints C with their corresponding Lagrange multipliers A. 

The principles of Rayleigh's dissipative function are described in some 
detail in textbooks on analytical dynamics, e.g. (l4|, |l5|| . Introduction of 
F = \likViVk-, with the velocities v = x and friction constants 7, al- 
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lows the unified description of Lagrangian systems where the total energy is 
dissipated at a rate 4r.E = —2F, via the variational equations 

A r A r 

— Ldt-— Fdt = 0. 
ox J ov J 

For instance, for a damped one-dimensional motion, 1Z = L — F = ^mx 2 — 
U(x) — and the variational equation correctly gives 

m x — 7 x + dU jdx = 0. 

This method has been proven to remain valid in problems involving an 



arbitrary set of Lagrangian constraints [16, 17| and promises a much greater 



mileage in non-equilibrium statistical mechanics of random networks. 



2 Tube model and equations of motion 

In order to describe the sliding motion of a polymer segment along the prim- 
itive path, we define the following variables, see Fig. 0(a). Let us consider a 
polymer chain of the arc length L. The arc distance from one end, s £ [0, L], 
is used to specify a position along the chain. The conformation of the poly- 
mer at time t is given, as a function of s, by the segment coordinate r(s,t). 
The chain is constrained within a tube of radius a - similarly, we take the 
conformation of the primitive path, at time t, given by the coordinate of 
tube axis R(n,t), where n £ [0,Z] is used to specify a position along the 
primitive path. Now, Za denotes the contour length of the primitive path 
and na is the arc distance from one end along the tube axis (a variable anal- 
ogous to s for the chain itself). An essential variable in the problem is the 
number of chain segments s contained within one step of the primitive path 
n. The relation between the position s on the polymer and the correspond- 
ing position n on the primitive path is then given by a function n = n(s,t). 
Note that r(s,t) and n(s,t) are dynamical variables of the polymer chain, 
while the tube conformation R(n, t) is given externally, by the surrounding 
constraints. 

In order to derive the over-damped equations of motion, we construct 
the Rayleighian 1Z as a functional of r, r, R and R, omitting the inertial 
effect of kinetic energy so that 1Z = —U — F. The potential energy term U, 
describing the Wiener random walk constrained around the path R(n, t) of 
a fixed length, is given by 

t>- z \- « 

Jo OS 

Here, b is the step length of the polymer chain and k^T is the Boltzmann 
temperature energy scale. The potential term containing q 2 represents the 



U 
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effect of confining the polymer segment r(s) around the primitive path 
through a harmonic potential. The last term in the eq. ([!]) arises from 
the constraint on the primitive path, n(L, t) — n(0, t) = Z, with a Lagrange 
multiplier A. 

The dissipative function F in the Rayleighian is determined by two rel- 
ative velocities - the movement of a chain segment with respect to its envi- 
ronment 

'dr dR 



dt dt 



s: fixed 



dR 

dn 



dR 

dt 



(2) 



and the slip motion of the chain segment equilibrium reference point along 
the primitive path 



dR 

~dt 



s: fixed 



dR 

~dt 



dR 

dn 



(3) 



where the local velocity of polymer segment and the rate of the primitive 
path evolution are given by, respectively, 



dr 

dt 



and 



dn 
~dt 



Two different friction constants are then introduced to describe the dissipa- 
tion corresponding to these two relative velocities, eqs. @ and (||), as they 
appear in the dissipative functional F: 



2 s 



ds v 



dR 

dn 



dR 

~dt 



1 



ds 



fdR 

V dn 
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V n ) ■ 



(4) 



The equations of motion are obtained from the Rayleighian 1Z through the 
variational equations 



6F SU 
5v 5r 







and 



SF 

5v n 



SU 
5n 



0. 



(5) 



The first of equations @ gives, after evaluating the functional derivatives, 



((—- ^i^L _ ^) + 3kBT 



V dt dn dt 



Of 



(6) 



The second varuiational equation (||) gives 



dR fdr dRdn dR\ fdR\ 2 dn 3k B T 2/ f dR\ 



which can be simplified to 

2 



\ dn ) dt 



dn J 



fdRy dn 3k B TdR 



\dn J dt 



dn 



d 2 r 

ds 2 



(7) 
(8) 
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with the help of eq. . 

A natural dynamical variable in our problem is not the absolute chain 
segment position r(s,t) but the position relative to the primitive path, the 
fact we have already exploited in eq. @ for the velocity. Accordingly, we 
introduce the transverse excursion 



f{s,t) = r(s,t) - R[n(s,t),t]. 



(9) 



The reason for using of the relative-position variable f(s,t) is to separate 
the fast and slow contributions to the dynamic equations. In particular, 
examine the eq. (||) for evolution of the primitive path: the equation of 
motion for n(s,t). If we neglect the effect of fast transverse fluctuations 
around the primitive path, f, this equation transforms to: 



fdR\ 2 dn 



t 3k B T dR 

\dn J dt b dn 



d 2 R 

ds 2 



0. 



(10) 



After one identifies the full derivatives 



dR 

~ds~ 



dRdn ^ d 2 R _ d 2 R fdn\ 2 dRd 2 n 
dn ds ds 2 dn 2 \ds J dn ds 2 



and makes the according substitution, the part of the dynamic equation (||) 
describing the slow evolution of the primitive path takes the form 



dn 3k B T 
~dt + bv 



dR d 2 R (dR 

dn dn 2 \ dn 



dnV 
~d~s) 



d 2 n 
ds 2 



0. 



:n) 



Note that, because s is measured in units of length, the dimensionality of 
v (and C) is such that the product (bv) is an ordinary friction constant 
measured, e.g., in J s/m 2 . In equilibrium, when no external deformation is 
applied to the system, the primitive path Rq is given by the static solution 
of the eq. (0) 



dRo d 2 Ro fdRoY 2 I 'dn 



dn dn 2 \ dn 



V ds 



d 2 n 
~ds^ 



0. 



(12) 



Since the equilibrium step length of the primitive path a is assumed to be 
fixed, one finds 



dRo 



dn 



dRo d 2 Ro 



0. 



dn dn 2 

With this orthogonality relation, the equilibrium condition (|l^) is simply 

d 2 n n , , Z 



(13) 



ds 2 



n (s 



L 



(14) 
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a homogeneous distribution satisfying the boundary conditions no(0) = 
and no(L) = Z. Let us note here, although a more extended discussion will 
follow later in the paper, that a limit of completely uncrosslinked system 
requires a non-trivial analysis. With respect to eq. (jlj) one simply needs 
to remember that Z and L are not independent parameters, in fact a direct 
proportionality Lb = Za 2 holds. Therefore, even for melt of uncrosslinked 
chains, or in the limit L — > oo, the constant ratio L/Z is always reflecting 
the average degree of entanglement, giving the arc distance between the 
constraints. 



3 The linear response 

We now consider the situation when an external deformation is suddenly 
applied to the network, and calculate the linear dynamic stress response. 
The principal feature of this analysis is the computation of non-equilibrium 
evolution of the tube surrounding the polymer chain. First, the primitive 
path trajectory, R(n,t), determined by local restrictions imposed by other 
chains in a dense system, is assumed to be deformed affinely with the whole 
body: 

R(n,t) =1 (t) -Ro(n) (15) 

with E (t) the deformation gradient tensor. For example, in a fixed coordi- 
nate system, E is given by 



1 Ai 2 
1 
1 



and 



Ai 
A 2 
A 3 



(16) 



for simple shear deformation and elongation, respectively. Since we shall be 
looking for the linear response, let us identify the small strain e (i) through 
E a /?(t) = $ a /3 + e a p(t). For the simple shear, only one component of strain 
tensor is present, e a p = eu a gp, where u is the displacement direction, g is 
that of a gradient and the amplitude e<l. For the extension, one obtains 
£q/3 = If the material is assumed fully incompressible, the imposed 

principle extension e = £3 -C 1 results in the other two diagonal elements 
equal to each other by rotational symmetry and £1 = £2 = 

It is necessary to separate the small non-equilibrium correction to the 
primitive path coordinate, n\, defined by n(s,t) = no(s) + m(s,t). Sub- 
stituting the affine deformation condition (^) and using the fixed-step and 
orthogonality relations (fllf), the linearised dynamic equation ( |TT|) takes the 
form of one-dimensional diffusion equation for a scalar function m(s,t): 



bv dn\ d 2 n\ 



3k B T dt ds 2 



(17) 
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with the external force 



ZydRo ,«-» j ~ d 2 R ( L \ d 



e 



aZ, 7 9n <9ri 2 \a 2 Z 1 ds 



dR 0a dR p 



-e 



ds ~ aP ~ ds 



(18) 

The diffusion equation of motion for the primitive path, eq. (17), is 
similar to what one frequently finds in the statistical theory of polymers, cf. 
|J. The diffusion constant of the homogeneous equation is D n = Sk^T/bv. 
In our case, since we shall be considering an oscillating imposed strain and 
looking for a response at a given frequency u, it is convenient to introduce 
the equilibrium Green's function that satisfies 



bv d 



2 



1W- 



G(s,s') = 5(s-s') (19) 



3k B T ds 2 ^ 
and is explicitly given by 

e -^ s - s '\ -= ( e - 7(W) sinh 7 s + e" 7S ' sinh 7 (L - s)\ 

simiTL I. > 



sinh 7L 

(20) 

with the characteristic parameter, of the inverse-length dimensionality, 



7 = J—- (21) 

This Green's function satisfies the required chain reversal symmetry and the 
boundary conditions 

G(s, s') = G(s', s); G(s, 0) = G(s, L) = 0. 



The dynamical solution of the eq. ( |17| ) with the perturbation ^(s, t) and the 
boundary conditions ni(0, t) = n\(L, t) = is then given by the convolution 

ni (s,t)= [ L ds'G(s,s')^(s',t) . (22) 
Jo 

The main dynamic equation (17), with the tensile force ^(s,t) due to effec- 
tive elongation of the primitive path under deformation, lead to a non-trivial 
relaxation behaviour - the key result of this paper. 

Once the linear correction to the primitive path is obtained [ni(s,t) 
generated by an imposed strain e (£)], one can compute the dynamic stress 
response. In a first approximation, the macroscopic stress is calculated 
through the correlation function along the chain (cf. 0), 

3/cbT f L /dr a dr \ 

where c x denotes the number density of polymer chains in the network (a 
quantity directly related to the crosslink density). In this equation, r(s) is 
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the absolute position of chain segment given by r(s, t) = R[n(s, t),t]+f(s, t). 
By neglecting the small fast-fluctuating variable f we obtain 

r(s) « Ro[n (s)}+ e (t) ■ Ro[n (s)} + ^f^ni(^), (24) 

up to the first order in e . Accordingly, one obtains three distinct contribu- 
tions to the average stress tensor, 



(7, 



a/3 



U a/3 + °a/3 + a a/3 ' 



corresponding to the undistorted primitive path Ro in the eq. (|24|), the 
equilibrium response to the imposed strain e and the relaxation part deter- 
mined by the sliding motion of polymer segments along the tube ni(s,t), 
respectively: 



a 



a 



a 



(0) 

a/3 

(1) 

a/3 

(2) 

a/3 



3k B T 
: b 

3k B T 
' b 

3k B T 



ii \ ds ds 

L 



ds 







^di?o,, d/?c, : . 



£/3 7 + 



di? /3 di?Qo 



ds ds / \ ds 

L ds / djR Q" djR o/3 5ra i( s ^) 
\ ds ds 



ds 



0S 



(25) 
,(26) 
(27) 



The first term is the small correction to the isotropic hydrostatic pressure 
in the system, determined by thermodynamic conditions (see, e.g. ||); 
we shall not discuss it any longer here. The second contribution, cr M, 
is the equilibrium entropic response of a chain confined to the primitive 
path random walk to an imposed strain; this is a time- independent effect. 
All leading-order effects of relaxation are concentrated in the third term, 
cr ( 2 ) , where the movement of chain segments is calculated with the help of 
equation (§3). 



4 The complex modulus G*(u) 

We now turn to particular cases of imposed small deformation and, first, 

consider the oscillating simple shear of the sample, cf. E s of eq. (|l6|); the 
corresponding small strain tensor is e a p = eu a gp e lw *. In order to evaluate 
averages in equations fl2q) and (|27]) one needs to recognize that, for the 
equilibrium primitive path Rq(jiq) one has, for instance, 

the pair average of the primitive path tangent vector t = (d-Ro/dn), which 
is assumed constant for all chain segments (labeled by s) belonging to the 
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same step n. In this way, the contour integral over the chain segments is 
divided into a sum over steps of the primitive path: 



ds 



z 

E 

n=l 



(n) 



4(n-l) 



ds. 



For the static stress a W this division does not matter because the integrand 
in (061) is not a function of s. One thus obtains 



a 



(i) 



3fc B T Z 2 1 2 . Nr 
Cx— £ ^-a £WM,3 + upg a \ 



2c x Zk B Te { ^, 



(29) 



where e w is the symmetric part of the imposed shear and the the relation 
Lb = Za 2 for the end to end distances of a chain and a primitive path is 
used. This equilibrium stress response is characterised by the effective rub- 
ber modulus n = 2c x Zk B T, proportional to the total number of constraints 
for the average chain. Note that because the chain of length L is considered 
to have (Z — 1) entanglement points and one cross link, the product c x Z 
can be written as c x + c e , involving the density of entanglements c e . 



Substituting the convolution relation (22) with the external force fy(s,t) 
given by the eq. (^) and integrating by parts to switch the s-derivative, we 
obtain for the relaxation part of stress: 



r (2) 
'a/3 



3k B T 



L 

aZ 



L dsds , d 2 G 
1 1 ds ds' 



dR 0a dRna dRo a dR 



10/3 



•0 b 



dn dn dn' dn' 



For simple shear, with vectors it-Lfif, only the four-term average gives a non- 
zero result, (t a tpt a t h ) = -^(^[Sapdab + 8 aa 8p b + 5 ab 5 ab ] at n = n'. One 
thus obtains, after replacing the integrals of full derivatives of the Green's 
function G(s, s'), 



a 



(2) 

a/3 



2k B T 

: 5b 



2k B T 



n=i 



#(n) 



Kn-1) 



(31) 



i) 



aZ 

T 



where 



and 




(S) 



(*)> 



5 - 2e- z? csch(Z0 sinh 2 (^) - - tanh(|jj32) 

2 Z 2 



(33) 



For a simple fixed- frequency oscillating strain input, the linear complex mod- 
ulus G*(u) for shear deformation is defined by the combination of the two 
relevant contributions to the stress tensor: 



G*(lv) = 2c x Zk B T 



1 



(34) 



Cab- (30) 
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As should be expected for a classical rubber, one obtains exactly the same 
expression for the linear modulus in response to the imposed extension e a p = 
£ a 5 a pe lujt . That is, an explicit calculation gives G*(u) = G*(u), where the 
latter is produced by the sum of two contributions to the normal stress, 



r (1) 

a/3 

r (2) 

T a(3 



2c x Zk B T£°J af3 $te e iujt 

2 [1 

--c x Zk B Te a 5 af 3?fte - 
5 U 



(35) 
(36) 



where c x Z 
£l = £2 = ~ 



= (c x + c e ) and the linear diagonal extensions are related by 
ie3 if the fully incompressible case is considered. 



5 Discussion 

In eq. (^J), the complex modulus G*(u) is given as a function of the com- 
plex non-dimensional parameter £, which is defined in eq. (p3|). Since this 
parameter contains the frequency to, the rest of it represents a time scale. 
This characteristic time, which we denote r e , is determined by the diffu- 
sion constant D n of the homogeneous equation of motion (|l7| ) for n\. For 
D n = 3k B T '/bv we obtain 



2D r , 



bv 
6k B T 



(37) 



which is the time for a density fluctuation n\ to diffuse a distance L/Z, a 
step length of the primitive path or the arc length of the polymer between 
the entanglements. This time scale is also found in the theory of polymer 
melts, cf. |^], where it separates the regimes of free Rouse dynamics and that 
constrained by the tube. The relation to the Rouse time is established when 
one recognises the product (bv) as the proper friction coefficient: r e ~ tr/Z 2 . 

The parameter £ can now be expressed as £ = (1 + i)^/T c uj, which also 
means that r e is one of the characteristic time scales of the complex modulus 
G*(uj). Examining the eq. d32| ) one can identify the second characteristic 
time appearing in the non-dimensional exponent Z^ ~ ^trlJ, with tr = 
L 2 /2D n the time for the slip- motion to diffuse all of the chain contour length 
between the crosslinks. 

The following limiting cases of non-dimensional function /(£, Z) can be 
easily seen from the definition (|32"|): 





Z -» 1 (C C /C X « 1) 


Z — > 00 


£^0 


^ — > 00 




^ L ) l+cosh(^) 


l-e-« 







which leads to a number of important limiting cases for the storage and loss 
parts of the complex modulus G*(uj) = G' + iG" : 



(arb. Z) 



G'(u) 



G"(u) 
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Figure 2: (a) Plots of storage modulus G'(lo), in absolute units, taking the 
value c x k B T = 10 5 J/m 3 (a typical rubber-elastic energy scale), against r e uj 
for Z =1,2,3,4 (increasing curves). The growth of the zero-frequency value 
Ho and a difference between it and the plateau value at r e uj 3> 1 are in 
evidence, (b) Plots of loss modulus G"(u>), for the same set of parameters 
and Zs, show an increasing effect of relaxation for more entangled network 
strands. 

The storage and loss moduli G' and G" are plotted against the reduced 
frequency in Figs. |2|(a,b), in absolute units: taking a typical for rubbers value 
c x k B T ~ 10 5 Pa. An important relation between the storage modulus value 
at zero frequency, the static modulus no, and that in the high-frequency limit 
(the "rubber plateau" value at T e uj 3> 1) is characterised by a universal ratio 

G"|t c u,»i 5Z 1 + (c e /c x ) 



G"|^o 4Z + 1 l + f(c e /e x ) 

Figure ||(a) makes this point more explicit by plotting the non-dimensional 
reduced storage modulus G' /2c x Z k B T, cf. eq. (|34|). At Z — > 1 there is 
no variation of G' with frequency. There is also no dissipation in this limit 
of high dilution, as shown by the plots of loss factor tan 5 in Fig. |3|(b). 
When entanglements dominate the rubber-elastic response, at Z > 1 (per- 
haps a more realistic case in typical rubber) one finds a 20% increase in the 
storage modulus for deformations on a time scale shorter than r e . There 
is a substantial dissipation effect in this limit of high entanglement, with 
the high-frequency tail of tan<5(o;) decaying rather slowly, according to the 
square-root asymptotic dependence, cf. eqs. (|38|). Note that, the characteris- 
tic time of network strand relaxation is rather short: in the reptation theory 
of melt dynamics one obtains the hierarchy of time scales r e <C tr < 
(the disengagement time, not present in an ideal crosslinked network) with 



tr/Z and ~ 3Zt b . Taking just a few recent cases [18[|, out a large 
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Figure 3: (a) Plots of reduced storage modulus G'/2c x Zk^T (non- 
dimensional) for Z =1,2,3,5 and 100 (decreasing curves), (b) Plots of the 
loss factor tan S = G"/G', for the same set of Zs (increasing curves). Note 
a slow decay of tan 5 (and increase of G') at t c uj 3> 1, which is due to the 
square-root approach to saturation values in eqs. (|38|). 



0.035 




0.02 0.04 0.06 0.08 0.1 

Reduced frequency T e C0 



Figure 4: Plots of the loss factor tan 5, for t c lj <C 1 and the increasing 
values of Z = 2, 4, 8 and 100. At moderate Z the linear regime is observed 
(the calculated slope is plotted as a dashed line for each Z). At large Z ^> 
(feu;) -1 / 4 one finds an apparent square-root scaling in a wide range of low 
frequencies. 



number of rheological studies of polymer melts, one obtains, very crudely, 
for polystyrene (PS) r e ~ 10~ 7 s and for poly(dimethyl) siloxane (PDMS) 
r e ~ 10~ 8 s. 

As a result, in some experiments, it may be difficult to access the high- 
frequency end of the predicted relaxation behaviour. It is, therefore, inter- 
esting to examine the asymptotic low-frequency behaviour of the loss factor 
for the cases of Z > 1 and Z> 1, when the number of entanglements per 
network strand becomes large. Equations ( |38|) indicate that at r e u> — ► 
the loss factor is linear: tan<5 ~ (Z — \)Z/(\2Z + 3)(r e o;). However, as Z 
increases so that Z S> (r e ti;) _1//4 , a different regime becomes apparent - see 
Fig. [|. In fact, this intermediate regime 1 S> t c uj ^> 1/Z 4 is best fit by an 
effective power-law tan 5 ~ 0.09(r e u;) ' 45 , nearly another square-root. 

The present calculation produces a static limit (u> — ► 0) which is exactly 
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equal to the result of Doi and Edwards M , derived by assuming that the de- 
formation of the primitive path is affine and that the tensile force is constant 
along the chain. Such a correspondence occurs because (i) the polymer den- 
sity along the tube is uniform for u = and (ii) in the presence of crosslinks 
at the both ends of the polymer there can neither be the contour length 
relaxation nor the disengagement relaxation. 

Ball, Doi, Edwards and Warner (BDEW) derived the formula for the 
elastic free energy within the sliplink model for entangled polymer networks, 
which produced the following expression of the static modulus: 



G' = 2k B T 



(39) 



where, in the notation of @, N c and N s denote the numbers of crosslinks 
and slip links, respectively. An additional parameter r\ characterises the 
amount of chain which can slip through a sliplink and, according to BDEW, 
is given by 

4 /(JV C + JV,)A* 2 



Here, L/ (N s + N c ) is the arc distance between two sliplinks, which is analo- 
gous to L/Z, the equilibrium step length of the primitive path, in the present 
paper, and A is the allowed slip arc length. By identifying that c x = N c 
and c e = N s for unit volume, we can see that the limits G'(r e u; ^> 1) and 
G'(r e uj <C 1) of the present calculation correspond to the following values 
of A in the BDEW result (^): At u — > oo one obtains A = (no slip- 
ping past the constraint) and at uj — > the slipping occurs to the amount 
A ~ 0.3(L/Z). 

One should note that a limiting case of zero crosslinking density, a corre- 
sponding entangled polymer melt, is not trivially achieved by setting c x — * 
in the consistent dynamical theory. At sufficiently short times, or frequency 
above t^uj ~ 1, the dynamic mechanical response is controlled by the chain 
constraints which are assumed to have the fixed density c e ~ Z/Ld 2 (for 
c x ~ 1/^chain = l/Ld 2 with d 2 the cross-section). This, and a higher- 
frequency regime when r e u > 1, is where the present model of sliding 
motion along the tube finds its application. The regime of chain disen- 
tanglement in a corresponding melt could be achieved, in a simple way, by 
setting phenomenologically an effectively decreasing entanglement density 
c| ff ~ exp (— t/Td), vanishing completely at t 3> due to reptation, see Q. 

To conclude, we have developed a first-approximation description of 
polymer dynamics in an ideal permanently crosslinked network under small 
deformations. This results in a general expression for linear complex modu- 
lus G*(u), which shows a characteristic single-relaxation time behaviour in 
its frequency variation. The low-frequency limits of storage and loss moduli 
have the classical dependence on frequency, square and linear, respectively. 
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We obtain a rather slow, square-root variation of G' and G" in the limit of 
high frequency, at t c uj 3> 1. The relaxation time scale r e , resulting from the 
description in terms of one-dimensional diffusion of chain segments along the 
distorted primitive path, is a quantity closely related to the characteristic 
entanglement (tube formation) time in a corresponding melt, although the 
topological sense (and, hence, the density of) entanglements changes when 
the chains are permanently crosslinked into a network. 

One of the principal achievements of this work is its consistent use 
of Rayleighian dissipative function approach. This method is naturally 
adopted for the treatment of fluctuations, relaxation and constraints, and 
could be useful for many other problems of polymer dynamics. 
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FIGURE CAPTIONS 



FIGURE 1. (a) The scheme of tube model and the relevant variables: the 
polymer conformation r(s, t) and the tube primitive path R(n, t), where the 
function n(s, t) specifies the position along the fixed-length path and s spec- 
ifies the position along the chain. In a network, the chain is permanently 
crosslinked at both ends of the tube. (b) The "expected" behaviour of 
the dynamic storage modulus G'(oj): For an uncrosslinked melt the static 
modulus is zero, but reaches a rubber plateau at a finite frequency due to en- 
tanglements. On increasing concentration of permanent crosslinks the static 
modulus ^'1^=0 = fJ-o increases; the plateau modulus has contributions from 
both crosslinks and entanglements. At a much higher frequency all curves 
should converge to the glass plateau, which is not the subject of this paper. 

FIGURE 2. (a) Plots of storage modulus G'(ll>), in absolute units, 
taking the value c x /cbT = 10 5 J/m 3 (a typical rubber-elastic energy scale), 
against t c uj for Z =1,2,3,4 (increasing curves). The growth of the zero- 
frequency value fio and a difference between it and the plateau value at 
T e uj ^> 1 are in evidence, (b) Plots of loss modulus G"(u), for the same 
set of parameters and Zs, show an increasing effect of relaxation for more 
entangled network strands. 

FIGURE 3. (a) Plots of reduced storage modulus G' /2c x Z ksT (non- 
dimensional) for Z =1,2,3,5 and 100 (decreasing curves), (b) Plots of the 
loss factor tan£ = G"/G', for the same set of Zs (increasing curves). Note 
a slow decay of tan S (and increase of G') at t c uj ^> 1, which is due to the 
square-root approach to saturation values in eqs. (|38[). 

FIGURE 4. Plots of the loss factor tan<5, for r e uj <C 1 and the 
increasing values of Z = 2, 4, 8 and 100. At moderate Z the linear regime 
is observed (the calculated slope is plotted as a dashed line for each Z). At 
large Z 3> (r c u;) _1 / 4 one finds an apparent square-root scaling in a wide 
range of low frequencies. 
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